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Abstract: In iKand (|2015f) we calculated the acceleration of cosmic-ray electrons and the ensuing radio 
synchrotron emission at weak spherical shocks that are expected to form in the outskirts of galaxy clusters. 
There we demonstrated that, at decelerating spherical shocks, the volume integrated spectra of both 
electrons and radiation deviate significantly from the test-particle power-laws predicted for constant planar 
shocks, because the shock compression ratio and the flux of injected electrons decrease in time. In this 
study, we consider spherical blast waves propagating into a constant density core surrounded by an 
isothermal halo with p oc in order to explore how the deceleration rate of the shock speed affects the 
radio emission from accelerated electrons. The surface brightness prohle and the volume-integrated radio 
spectrum of the model shocks are calculated by assuming a ribbon-like shock surface on a spherical shell 
and the associated downstream region of relativistic electrons. If the postshock magnetic field strength is 
about 7 /rG, at the shock age of 50 Myr, the volume-integrated radio spectrum steepens gradually with 
the spectral index from ainj to a-mj + 0.5 over 0.1-10 GHz, where ainj is the injection index at the shock 
position expected from the diffusive shock acceleration theory. Such gradual steepening could explain the 
cu rved radio spec t rum o f the radio relic in cluster A2266, which was interpreted as a broken power-law 
by iTrasatti et ahl (l2014ll . if the relic shock is young enough so that the break frequency falls in around 
1 GHz. 
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1. INTRODUCTION 


Cosmological hydrodynamic simulations have shown 
that shock waves may form due to supersonic flow 
motions in the baryonic intracluster medium (ICM) 
during the formati on of the large scale structure in 
the Universe (e-K - - iR-vu et al.1 2003t Kang et H] l2007t 
IVazza et al.ll200il ISkillman et al. 201 ill . The time evo¬ 
lution of the spatial distri bution of such sho cks in nu¬ 
merical simulations fe.g.. I^zza et al1l2012f) indicates 
that shock surfaces behave like spherical bubbles blow¬ 
ing out from the cluster center, when major episodes of 
mergers or infalls from adjacent filaments occur Shock 
surfaces seem to last only for a fraction of a dynamical 
time scale of clusters, i.e., O.ltdyn 100 Myr. This 
implies that some cosmological shock waves are associ¬ 
ated with merger induced outflows, and that spherical 
shock bubbles are likely to expand into the cluster out¬ 
skirts with a decreasing density profile. 

Observational evidence for shock bubbles can be 
found at the so-called “radio relics” detected in the 
outskirts of galaxy cluster, which are interpreted as 
synchrotron emitting structures containing relativistic 
electr ons accelerated at weak I CM shocks (M„ ^ 2 — 4 


e.g.. Ivan Weeren et al.1 2^^ 20111: Nuza et 1201 
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iFeretti et al.ll201 Brunetti fc Jone^boidll . Some radio 

relics, for instance, the “Sausage relic” in galaxy clus¬ 
ter GIZA J2242.8-1-5301 and the “Toothbrush” relic in 
galaxy cluster IRXS J0603.3-1-4214, have thin arc-like 
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shapes of ^ 50 kpc in width an d ^1 — 2 Mpc in length 
(|vaii Weeren et al1l2010L l2012ll . They could be repre¬ 
sented by a ribbon-like structure o n a spherical shell 
proje cted onto t he sky plane (e.g., Ivan Weeren et al.l 
l2010HKang et al.ll2012ll . 

As shown in Fig. 1, the viewing depth of a relic shock 
structure can be parameterized by the extension angle 
2'(/i, while the physical width of the postshock shell of ra¬ 
diating electrons is mainly determined by the advection 
length, AGdv(7e) ~ U 2 ■ min[tage,trad(7e)], where U 2 is 
the flow speed behind the shock and trad(7e) is the ra¬ 
diative cooling time scale for electrons with Lorentz fac¬ 
tor, 7 e. It remains largely unknown how such a ribbon¬ 
like structure can be formed by the spherical outflows 
in galaxy clusters. We note, however, a recent study by 
IShimwell et al.l (|2015ll who suggested that the uniform 
arc-like shape of some radio relics may trace the un¬ 
derlying region of pre-existing seed electrons remaining 
from an old radio lobe. 

According to diffusive shock acceleration (DSA) the¬ 
ory, cosmic-ray (GR) particles can be generate d via 
Ferni i 1st order process at collisionless shocks (jBelll 
119781: [Prurvl Il983t iMalkov fc Drur^ l200lll . The test- 
particle DSA theory predicts that the GR energy spec¬ 
trum at the shock position has a power-law energy spec¬ 
trum, N{E) oc where s = (ct -|- 2)/(ct — 1) and 

a = P 2 /Pi is the shock compression ratio. Hereafter, 
we use the subscripts ‘I’ and ‘2’ to denote the condi¬ 
tions upstream and downstream of shock, respectively. 
Then the synchrotron radiation spectrum due to CR 
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Figure 1. Schematic diagram showing a spherical shock surface (red) and the postshock electron distribution behind the shock 
adopted in this study. The shock surface is modeled as a Mpc-long ribbon on the spherical surface whose viewing depth 
can be defined by the extension angle 2ip. The depth of the postshock electron distribution behind the shock is determined 
by the advection length, Aiadv ~ U 2 • tage, for low energy electrons or the cooling length, Aicooi(7e) ~ U 2 ■ tradije) for high 
energy electrons. The radio emissivity from CR electrons, jv{d) is calculated using the DSA simulation results in Paper I 
(where d is the distance from the shock). The surface brightness is calculated by integrated jv{d) over the path along the 
line of sight (see Eq. my If viewed from th e top as indicated here, the projected image will have an arc-like shape similar 
to the Sausage relic in CIZA J2242.8-1-5301 llvan Weeren et ^l20inh . 


electrons at the shock location has a power-law form of 
jv{xs) oc where ainj = (s — l)/2 is the injection 

index. Moreover, the volume-integrated synchrotron 
spectrum downstream of a planar shock becomes a sim¬ 
ple power-law of cx with A^, = ainj + 0.5 

above a break frequency, since electrons cool via syn¬ 
chrotron an d inverse-Co mpton (iC) losses behind the 
shock (e.g., iKand I20T1II . Such spectral characteris¬ 
tics are commonly used to infer the shock properties 
(i.e., Mach number) of observed radio reli c shocks (e.g., 
Ivan Weeren et al.llioiOl : IStroe et 701201411 

In iKand (|2015D (Paper I), we calculated the electron 
acceleration at spherical shocks similar to Sedov-Taylor 
blast waves with Mg 2.5 — 4.5, which expand into a 
hot uniform ICM. We found that the electron energy 
spectrum at the shock location has reached the steady 
state defined by the instantaneous shock parameters. 
So the spatially resolved, synchrotron radiation spectra 
at the shock could be described properly by the test- 
particle DSA predictions for planar shocks. However, 
the volume integrated spectra of both electrons and ra¬ 
diation evolve differently from those of planar shocks 
and exhibit some nonlinear signatures, depending on 
the time-dependent evolution of the shock parameters. 
For instance, the shock compression ratio a and the 
injection flux of CR electrons decrease in time, as the 
spherical shock expands and slows down, resulting in 
some curvatures in both electron and radiation spec¬ 
tra. 


Magnetic fields play key roles in DSA at collisionless 
shocks and control the synchrotron cooling and emis¬ 
sion of relativistic electrons. Observed magnetic field 
strength is found to decrease from ~ 1 — 10 in 
the core region to "^ 0-1 ~ 1 /rC in the periphery of 
clusters (jFeretti et al.l l2012h . On the other hand, it 
is well established that magnetic fields can be ampli¬ 
fied via resonant and non-resonant instabilities induced 
by CR protons streaming upstrearn of st rong shocks 
(iBelil 119781: [Lucek fc Belli 120001: lB^[2n?)l . Recently, 
ICanrioli fc Soitkovskvl (l2014ll has shown that the mag¬ 
netic field amplification (MFA) factor scales with the 
Alfvenic Mach number, Ma, and the CR proton ac¬ 
celeration efficiency as {6B/B)‘^ ~ 3Ma{Pci, 2/Piu^). 
Here SB is the turbulent magnetic fields perpendicu¬ 
lar to the mean background magnetic fields, pi is the 
upstream gas density, and Per,2 is the downstream CR 
pressure. For t ypical cluster sho cks with 3 < Mg < 5 
and Ma '^10 (iRvu et al.l 120031) . the MFA factor due 
to the streaming stabilities is expected be rather small 
but not negligible, {SB/B)"^ ^ 0.3 — 1.5. However, 
it has not yet been fully understood how magnetic 
fields may be amplified in both upstream and down¬ 
stream of a weak shock in high beta ICM plasmas with 
Pp = Pg/PB 100 . 


Therefore, in Paper I, we considered several mod¬ 
els with decaying postshock magnetic fields, in which 
the downstream magnetic field, Bd{r), decreases be¬ 
hind the shock with a scale height of 100-150 kpc. We 
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found that the impacts of different Bd{r) profiles on 
the spatial distribution of the electron energy spectrum, 
Ne{r,^e), is not substantial, because iC scattering off 
cosmic background photons provides the baseline cool¬ 
ing rate. Any variations in iVe(r, je) is smoothed out in 
the spatial distribution of the synchrotron emissivity, 
jv^r), because electrons in a broad range of 7 e con¬ 
tribute to jv at a given frequency. We note, however, 
that the dependence of the synchrotron emissiv¬ 
ity {jv ot Nf.B^) can become significant in some cases. 
Moreover, any nonlinear features due to the spatial vari¬ 
ations of Nf.{r,je) and Bd{r) are mostly averaged out, 
leaving only subtle signatures in the volume integrated 
spectrum, J^. 

For the case with a constant background density (e.g., 
MFl-3 model in Paper 1), the shock speed decreases 
approximately as Ug oc In fact, it has not been 

examined, through cosmological hydrodynamic simula¬ 
tions, whether these shock bubbles would accelerate or 
decelerate as they propagate into the cluster outskirts 
with a decreasing density profile. In this study, we 
have performed additional DSA simulations, in which 
the initial Sedov blast wave propagates into the back¬ 
ground medium with pu oc r"”, where r is the ra¬ 
dial distance from the cluster center and n = 2 — 4. 
This effectively mimics a blast wave that expands into 
a constant-density core surrounded by an isothermal 
halo with a decreasing pressure profile. In these new 
runs, the spherical shock decelerates much slowly, com¬ 
pared to Ug oc and the nonlinear effects due to 

the deceleration of the shock speed are expected to be 
reduced from what we observed in the uniform density 
models considered in Paper I. 

Here we also have calculated the projected surface 
brightness, which depend on the three dimensional 
structure of the shock surfaces and the viewing direc¬ 
tion. Because of the curvature in the model shock struc¬ 
ture, synchrotron emissions from downstream electrons 
with different ages contribute to the surface bright¬ 
ness along a given line-of-sight (i.e., projection effects). 
So the observed spatial profile of Iu{R) is calculated 
by assuming the geometrical configuration described 
in Fig. I. Note that this model with a ribbon-like 
shock surface gave rise to radio flux profiles, Su{R), 
that were consistent with those of the the Sausage 
relic in CIZA J2242.8-1-5301 and the double relics in 
ZwC10008.8T521 5 (e.g., Ivan Weeren et abll2mni 120111: 
iKang et al.llMl . 

In paper I, we demonstrated that the spectral index 
of the volume-integrated spectral index increases grad¬ 
ually from Ai, = Oinj to = Oinj -1-0.5 over a broad 
frequency range, ^ (0.1 — 10)ubr, where the break fre¬ 
quency is Ubr ~ 0.5 GHz at the shock age of about 
50 Myr for the postshock magnetic fields ^ 7 pG (see 
Eq. 0). Here we have explored whether such a transi¬ 
tion can explain the broken power-law spectra ob served 
in the radio relic in A2256 (jTrasatti et al.lU014l1 . 

In the next section we describe the numerical cal¬ 
culations. The DSA simulation results of blast wave 
models with different postshock magnetic field profiles 


and with different background density profiles will be 
discussed in Section 3. A brief summary will be given 
in Section 4. 

2. Numerical Calculations 

In order to calculate DSA of CR electrons at spherical 
shocks, we have solved the time-dependent diffusion- 
convection equation for the pitch-angle-averaged phase 
space distribution function for CR electrons, /e(r’, p, t), 
in the one-dim ensional spherically symmetric geometry 
(ISkillinBi lI975^ . The test-particle version of CRASH 
(Cosmic-Ray Amr SHock) code in a comoving spherical 
grid was used (iKang fc ,IoneEll2?)?)^ . The details of the 
simulation set-up can be found in Paper I. 

For the initial conditions, we adopt a Sedov-Taylor 
similarity solution propagating into a uniform ICM 
with the following parameters: the ICM density, nup = 
10“^ cm“^, the ICM temperature, Ti = 5 x lO^K, 
the initial shock radius, Vgp = 0.78Mpc, and the ini¬ 
tial shock speed, Ugp = 4.5 x 10^ km s“^ with the 
sonic Mach number, Mgp = 4.3 at the onset of the 
simulation. The shock parameters change in time as 
the spherical shock expands out, depending on the up¬ 
stream conditions in the cluster outskirts. For iC cool¬ 
ing, a redshift z = 0.2 is chosen as a reference epoch, 
so Brad = 3.24 pG • (I -k zf = 4.7 pG. 

The synchrotron emissivity, jui'i'), at each shell is 
calculated (in units of erg cm“^ s“^ Hz“^ str“^), us¬ 
ing the electron distribution function, /e (r, p, t) , and 
the magnetic field profile, B{r,t), from the DSA sim¬ 
ulations. Then the radio intensity or surface bright¬ 
ness, Ii, (erg cm“^ s“^ Hz“^ str“^), is calculated by 
integrating j^, along the path length, h, as shown in 
Fig. 1: 

rhmax. 

R{R) = 2 / Md)dh. (I) 

Jo 

Here R is the distance behind the projected shock edge 
in the plane of the sky and d = Vg—r is the distance of a 
shell behind the shock, where (ug — = (ug — i?)^ -|- h^. 

The extension angle '0 = 10° is assumed in this study. 

The volume integrated emissivity, = f j^(r)dV, is 
calculated by integrating over the downstream vol¬ 
ume defined in Fig. I. 

The spectral indexes of jv{d), J^, and R{R) are de¬ 
fined as follows: 




A 


B^i-i^i+iiR) = 


d in j„{d) 
dinv ’ 
din Ji, 
dlnu ’ 
d\nI^{R) 
dinv ’ 


( 2 ) 

( 3 ) 

( 4 ) 


estimated between Ui and We chose the follow¬ 

ing three frequencies at the source, vi = 240 MHz, 
V 2 = 600 MHz, z /3 = 1.4 GHz in Figs. 2 and 4. Then 
the redshifted frequency for objects at a redshift z is 

I'obs.i = t'i/)! + z). 
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d or R (kpc) d or R (kpc) d or R (kpc) 

Figure 2. Spherical shock models with three different magnetic field profiles: MFl, MF2, and MF3. The results are shown 
at tage = 47 Myr, when rs = 0.96Mpc, Us = 3.4 x 10® km s“® and Mb = 3.2. In the upper two rows, spatial distributions of 
the synchrotron emissivity, Vijviir), and its spectral indexes, are plotted as a function of the postshock distance 

d from the shock surface. In the lower two rows, the intensity Uilv^iR) and its spectral indexes, are plotted 

as a function of the projected distance R from the shock. The frequency is Oi = 240 MHz (black solid line), 600 MHz (red 
dotted), 1.4 GHz (blue dashed), and 3 GHz for i = 1, 2, 3, and 4. The downstream volume of radio-emitting electrons is 
assumed to have the same shape as the one illustrated in Fig. 1 with the extension angle ip = 10°. Note that j^, and R are 
plotted in arbitrary units. 


Table 1. Model Parameters 


Model 

Pv.{r) 

BRr) 

MF1°^ 

pu ~ pO 

Bd = 7pG 

MF2 

pu = PO 

Bd OC p(r) 

MF3 

pu = PO 

Bd OC ^P{r) 

BDl*" 

Pu = PO 

Bd = 7pG 

BD2 

Pu OC 

Bd = 7pG 

BD2b 

Pu OC r~^ 

Bd OC ^P(r) 

BD3 

Pu OC r“'* 

Bd = 7pG 

BD3b 

Pu OC r“'* 

Bd OC ^JP{r) 


In fact MFl and BDl models are identical. 


3. DSA SIMULATION RESULTS 

3.1. Shocks with Different Magnetic Field Profiles 

As in the previous study of iKang et al.l (|2ni2ll . we adopt 
the postshock magnetic field strength ^2 ~ 7 fj,G in 
order to explain the observed width of the Sausage 
relic (AZ ~ 50 kpc), while the preshock magnetic field 
strength is chosen to be Bi ~ 2 — 3 fiG. Since Bi can¬ 
not be constrained directly from observations, it is ad¬ 
justed so that i ?2 becomes about 7 fiG after considering 
MFA or compression of the perpendicular components 
of magnetic fields across the shock. 

We consider several models whose characteristics are 
summarized in Table 1. As described in Paper I, we 
adopt the following downstream magnetic field profile, 
Bdir) in MF models, in which the upstream density 


























Radio Emission from Spherical Cluster Shocks 


5 



Figures. Evolution of spherical shock models with different background density profile, BDl (black solid lines), BD2 (red 
dotted), and BD3 (blue dashed): shock radius, Rs, shock speed, Ua, sonic Mach number, Ma, and the DSA spectral index 
at the shock, ainj. 


and temperature are uniform: 

MFl: Bi = 2 fiG k B2 = 7 fiG. 

MF2: = 3 ^iG, B 2 = Biy/TJf+2^^, 

& Bd{r) = B 2 ■ (F(r)/P 2 )^^^ for r < Vg. 

MF3: = 3 /rG, B 2 = Siy/l/3 + 2cr2/3, 

& Bd{r) = B2 ■ {p{r)/p2) for r < 

In Fig. 2 we compare these models at the shock 
age of 47 Myr. Note the evolution of the shock (i.e., 
rs{t) oc and Us{t) oc is identical in the three 

models. Synchrotron emissivity scales with the elec¬ 
tron energy spectrum and the magnetic field strength 
as jv oc Ne{'-ff.)B^. Of course, the evolution of Ne{^e) in 
each MF model also depends on the assumed profile of 
Bd{r) through DSA and synchrotron cooling. In MFl 
model where Bd{r) = B2, the spatial distribution of 
ju (d) at 240 MHz (black solid line) reveals dramatically 
the shock deceleration effects. It increases downstream 
behind the shock because of higher shock compression 
and higher electron injection flux at earlier time. On 
the other hand, j^{d) at 1.4 GHz (blue dashed line) de¬ 
creases downstream behind the shock due to the fast 
synchrotron/iC cooling of high energy electrons. 

The downstream increase of ji^{d) at low frequen¬ 
cies is softened in the models with decaying postshock 
magnetic fields because of dependence of the syn¬ 
chrotron emissivity. In MF3 model, in which the down¬ 
stream magnetic field decreases with the gas density 
behind the shocks, the spatial distributions of ji^{d) 
at all three frequencies decrease downstream. Thus 
the downstream distribution of ju{d) depends on the 
shock speed evolution, Us{t) and the postshock Bd{r) as 


well as the chosen frequency. In Paper I, we indicated, 
rather hastily, that signatures imprinted on synchrotron 
emission, ji/(r), and its volume integrated spectrum, J^, 
due to different postshock magnetic field profiles could 
be too subtle to detect. The main reason for this dis¬ 
crepancy is that we examined log j^{r) in Paper I (see 
Fig. 6 there), while vj^id) is plotted in Fig. 2. 

As can be seen in the lower two rows of panels in 
Fig. 2, the surface brightness profiles, Iv{R), are af¬ 
fected by the projection effects as well as the evolution 
of Us{f) and the spatial variation of Bd{r). For instance, 
the gradual increase of R just behind the shock up to 
i?inf <15 kpc is due to the increase of the path length, 
and its inflection point, Rmi = Ts(1 — cosip), depends 
on the value of ip. Beyond the inflection point, the path 
length decreases but ji, (d) may increase or decrease de¬ 
pending on Us(t) and Bd(r), resulting in a wide range 
of spatial profiles of Iv{R). 

Fig. 2 also demonstrate that the spectral indexes, 
and (i?), at all three frequencies 

decrease behind the shock and do not show significant 
variations among the different MF models other than 
faster steepening of both indexes for a more rapid decay 
of Bd(r). 

In summary, at Sedov-Taylor type spherical shocks 
decelerating with Ug oc the distribution function, 

Ne(r, 7 e), of low energy electrons increases downstream 
behind the shock. As a result, the spatial distributions 
of the radio emissivity, jv{d), and the surface bright¬ 
ness, Iu{R), at low radio frequencies (< 1 GHz) could 
depend significantly on the postshock magnetic field 
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R (kpc) R (kpc) R (kpc) 

Figure 4. Spherical shock models with different background density and magnetic field profiles: BDl (MFl), BD2, BD2b, 
and BD3b (from top to bottom panels). See Table 1 for the model parameters. Spatial distributions of the intensity vlv 
are shown at three shock ages, = 13 Myr, t 2 = 33 Myr, and ts = 67 Myr. The frequency is = 240 MHz (black solid 
line), U 2 = 600 MHz (red dotted), and ui = 1.4 GHz (blue dashed). The downstream volume of radio-emitting electrons is 
assumed to have the same shape as the one illustrated in Fig. 1 with the extension angle ip = 10°- The intensity, uli, ■ X, is 
plotted in an arbitrary unit, where the numerical factor, 1 < X < 5, is adopted in order to plot the quantities in the linear 
scale. 


profile. At high radio frequencies, such dependence be¬ 
comes relatively weaker, because the postshock width 
of high energy electrons is much narrower and so the 
magnetic field profile far downstream has less influence 
on synchrotron emission. On the other hand, the spec¬ 
tral indexes, and 5,^, are relatively insensitive to 
those variations. 

3.2. Shocks with Different Background Density Profile 

In BD models, we assume that the initial blast wave 
propagates into an isothermal halo with a different den¬ 
sity profile (r > Tg^i): 

BDl: puir) = po. 

BD2: pu{r) = po ■ 

BD3: pu{r) = po ■ {r/rg^i)~'^. 

Again the upstream temperature (Ti = 5 x 10^ K) 
is uniform, and Bi = 2 pG & i ?2 = 7 pG for these 
models. In fact MFl and BDl models are identical. 


We also consider BD2b and BD3b models, in which 
the downstream magnetic field profile is the same as 
MF2, i.e., Bd{r) cx y^P(r). 

In the so-call beta model for the gas distribution for 
isotherm al ICM, p(r) o c in the outskirts of galaxy 
clusters ( SarazinI Il988ll . So BD2 model corresponds 
to the beta model with pi r; 2/3, which is consistent 
with typical X-ray brightness profile of observed X-ray 
clusters. Recall that in our simulations the spherical 
blast wave into a uniform ICM is adopted for the initial 
conditions. So we are effectively considering a spherical 
blast wave that propagates into a constant-density core 
(i.e., p = Po ioT r < Vg^i) surrounded by a isothermal 
halo with p oc (for r > Tg^i). 

With the different background density (or gas pres¬ 
sure) profile, the shock speed evolves differently as it 
expands outward. Fig. 3 shows how the shock ra¬ 
dius, rs(t), shock speed, Us{t), the sonic Mach number. 
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Figures. Same models shown as in Fig. 4. Volume integrated emissivity, J^, and its spectral index, A^, are shown at three 
different shock age, ti — 13 Myr (black solid lines), t 2 = 33 Myr (red dotted lines), and is = 67 Myr (blue dashed lines). 


Ms{t), and the DSA spectral index at the shock posi¬ 
tion, Oinj = 3/[2 ((t — 1)], varies in time in BD models. 
As expected, the shock decelerates much slowly if the 
background pressure declines outward as in BD2 and 
BD3 models. In fact, the shock speed is almost con¬ 
stant in BD3 model. Hence these BD models allow 
us to explore the dependence of radio spectral proper¬ 
ties for a wide range of the time evolution of the shock 
speed. 

Fig. 4 compares the surface brightness profile, Iv{R), 
at vi = 240 MHz (black solid lines), V 2 = 600 MHz (red 
dotted), and U 3 = 1.4 GHz (blue dashed) in the models 
with different pu and B^. In BDl model, the shock 
parameters are: Rg = 0.83Mpc, Ug = 4.1 x 10^ km s“^ 
and Mg = 3.9 at ti = 13 Myr; Rg = 0.91Mpc, Ug = 
3.6 X 10^ km s“^ and Mg = 3.4 at ^2 = 33 Myr; Rg = 
l.OMpc, Ug = 3.1 X 10^ km s“^ and Mg = 2.9 at = 
67 Myr. Again, we can see that the surface brightness 
profile of I^{R) depends on the evolution of Ug{t) and 
the spatial profile of Bd{r). 

The width of the radio structure, Aladv ~ Gge^s/cr, 
increases with the shock age, while the amplitude of 


the surface brightness decreases in time (from left to 
right in Fig. 4). For high frequencies (> 1 GHz), 
however, the width asymptotes to the cooling length, 
Alcooi ~ irad( 7 e)'Ws/cr. So the postshock magnetic field 
strength can be obtained through trad( 7 e) from high 
frequency observations. At low frequencies (< 1 GHz) 
the surface brightness is affected by the evolution of 
Ug{t) and the spatial variation of Bd(r) as well as the 
projection effects, resulting in a wide range of spatial 
profiles of Iv{R)- This suggests that, if the radio surface 
brightness can be spatially resolved at several radio fre¬ 
quencies over ~ (0.1 — 10) GHz, we may extract the time 
evolution of Ug (t) as well as the shock age from low fre¬ 
quency observations. So it would be useful to compare 
multi-frequency radio observations with the intensity 
modeling that accounts for both shock evolution and 
projection effects. 

The comparison of R of MFl (BDl) and MF2 
models (see Fig. 2) demonstrates that the decaying 
postshock magnetic field could lessen the shock decel¬ 
eration signatures in the surface brightness profile. If 
we compare Ii, of BD2 and BD2b models in Fig. 4, 
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on the other hand, the difference in their profiles is 
much smaller than that between MFl and MF2 mod¬ 
els. This is because the the shock deceleration effects 
and the ensuing downstream increase of Nf,{r,"fe) are 
relatively milder in BD2 and BD2b models. 

Next Fig. 5 compares the volume integrated spec¬ 
trum, J[^, and its spectral index, at three different 
shock ages, tage = 13 Myr (black solid), 33 Myr (red 
dotted), and 67 Myr (blue dashed) for the same set of 
models shown in Fig. 4. 

For the test-particle power-law at planar shocks, the 
radio index Ai, is the same as ainj for v < ^'br, where 


r'br 


0.44GHz 


(sOMyr 


t’age 


-2 


B. 


e,2 


8.4 /iG 


-4 



(5) 

at the source. For v > z/br, A^ is expected to steepen 
to Ginj -1-0.5 due to synchrotron/iG cooling of electrons. 
Note that the smallest possible index is ainj = 0.5 at 
strong shocks, so A^, = ainj -1-0.5 > 1.0 above the break 
frequency. Also note that the observed (redshifted) 
break frequency corresponds to ;/br,obs = r'br/(l + z). 
If the shock age is Gge ~ 67 Myr and H 2 ~ 7 fj,G, the 
observed break frequency becomes about 120 MHz for 
objects at z = 0.2. 

As can be seen in Fig. 5, however, the transition 
from Ginj to Ginj + 0.5 occurs gradually over about two 
orders of magnitude in frequency range. If the shock 
is 30-70 Myr old for the shock parameters considered 
here, the volume-integrated radio spectrum is expected 
to steepen gradually from 100 MHz to 10 GHz, instead 
of a sharp broken power-law with the break frequency 
at t'br.obs- Such gradual steepening may explain why 
the volume-integrated radio spectrum of some relics was 
interpreted as a broken power-law with A^, < 1.0 at low 
frequencies. 

According to a recent observation of the relic in 
A2256, for example, the observed spectral index is 
Hggf® Rs 0.85 between 351 and 1369 MHz and in- 
creases to Aj°Ro° 1.0 between 1369 and 10450 MHz 
(|Trasatti et al.l [201411 . We estimated similar indexes 
for our Jy between 355 and 1413 MHz, and between 
1413 MHz and 10 GHz, using the DSA simulation re¬ 
sults. In BD2b model, for example, at tage = 33 Myr 
(red dotted line in Fig. 5), these spectral indexes are 
_ 0.85 and «:! 1.07. In BD3b model, 

again at Gge = 33 Myr, they are Ri 0.83 and 

^ 1413 ° ~ 1-03. So if the shock age or the electron ac¬ 
celeration duration is about 30 Myr, relatively young 
compared to the dynamical time scale of typical clus¬ 
ters, the curved radio spectrum around 1 GHz can be 
explained by DSA at cluster shocks. 

In the case of BD3b model, the spectral index be¬ 
haves very similarly to that of a plane shock case (see 
Fig. 4 of Paper 1), since the shock speed is more or 
less constant in time. Departures from the predictions 
for the test-particle planar shock are the most severe in 
BDl (MFl) model, while it becomes relatively milder 
in BD2 and BD3 models with decreasing halo density. 

As shown in Paper I, any variations in the spatial 


distributions of /e(r) and Bd(r) are smoothed in the 
volume integrated quantities such as J^. So signatures 
imprinted on the volume-integrated emission due to dif¬ 
ferent Bd{r) (e.g., between BD2 and BD2b models) 
would be too subtle to detect. 


4. SUMMARY 


We have performed time-dependent DSA simulations 
for cosmic-ray (GR) electrons at decelerating spheri¬ 
cal shocks with parameters relevant for weak cluster 
shocks: Ug ~ (3.0 — 4.5) x 10^ km s“^ and Mg ~ 
3.0 — 4.3. Several models with different postshock mag¬ 
netic field profiles (MFl-3) and different upstream gas 
density profiles (BDl-3) were considered as summa¬ 
rized in Table 1. Using the synchrotron emissivity, 
calculated from the CR electron energy spectra 
at these model shocks, the radio surface brightness pro¬ 
file, Iv(R), and the volume integrated spectrum, J^, 
were estimated by assuming a ribbon-like shock struc¬ 
ture described in Fig. 1. 

At low frequencies (< 1 GHz) the surface brightness 
is affected by the evolution of Us{t) and the spatial vari¬ 
ation of Bd{r) as well as the projection effects , result¬ 
ing in a wide range of spatial profiles of /^(i?) (see Figs. 
2 and 4). At high frequencies (> 1 GHz), such depen¬ 
dences become relatively weaker, because the postshock 
width of high energy electrons is much narrower and so 
the magnetic field profile far downstream has less influ¬ 
ence on synchrotron emission. 

For low frequency observations, the width of radio 
relics increases with the shock age as AZadv ^ageWs/cr, 
while it asymptotes to the cooling length, AUooi ~ 
Uad( 7 e)Ms/o', for high frequencies. If the surface bright¬ 
ness can be spatially resolved at multi-frequency obser¬ 
vations over ~ (0.1 — 10) GHz, we may extract sig¬ 
nificant information about the time evolution of Ug{t)^ 
the shock age, tage, and the postshock magnetic field 
strength, Bd{r), through the detail modeling of DSA 
and projection effects. On the contrary, the spectral 
index of R^R) behaves rather similarly in all the mod¬ 
els considered here. 

If the postshock magnetic field strength is about 
7 /rG, at the shock age of ^ 30 Myr, the volume- 
integrated radio spectrum has a break frequency, i/br 
1 GHz, and steepens gradually with the spectral index 
from ttinj to Ginj 4-0.5 over the frequency range of 0.1- 
10 GHz (see Fig. 5). Thus we suggest that such a 
curved spectrum could expla in the observed spec trum 
of the relic in cluster A2256 ( Trasatti et al.l[^014^ . 
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